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ABSTRACT 


Landsat-4 Thematic Mapper (TM) images recorded on computer 
compatible tapes (CCTs) are currently available from the SCROUNGE 
image processing system of the Landsat-4 Assessment System (LAS) 
project of the National Aeronautics and Space Administration 
(NASA) at Goddard Space Flight Center (GSFC) in the following 
tape formats: 1) raw rectified data (CCT-BT); 2) calibrated 

data (CCT-AT); and 3) geometrically corrected data (CCT-PT). 
These formats represent different steps in the process of 
producing fully-corrected TM data. CCT-BT images are rectified 
from telemetry format to image format, but are uncorrected radio- 
metrically and geometrically; they are generally used for 
internal transportation of data from one ground processing system 
to another, CCT-AT images have had data from faulty channels 
replaced and all data radiometrically calibrated to produce an 
archive image; they are available to researchers for radiometric 
characterization. The final products, the CCT-PT images, have 
been resampled by cubic convolution procedures to provide a 
geometrically corrected image using satellite ephemeris data. 

The CCT-PT image is the one to which all of the various radio- 
metric and geometric corrections have been applied; this is the 
product that is available to all users. 

Image files on the three tapes of the CCT-PT set are blocked 
four image lines per tape record, yielding a blocksize of 28672 . 
bytes. The CCT-AT images also has a blocking factor of four but 
a blocksize of 26624 bytes; the equivalent values for the CCT-BT 
images are five and 32000 bytes. 

Forward and reverse scans on all products are properly 
oriented east-west, but the CCT-BT and CCT-AT images have an 
alignment offset between scans. Original band-6 120m IFOV pixels 
have been replicated to four 30m IFOV pixels along a line in 
CCT-BT images, but the lines are not replicated to yield a full 
4x4 block of 30m IFOV pixels that are array compatible with the 
other bands until the CCT-AT image is produced. 

Saturated or very bright isolated pixels occurring in all 
bands result from apparent specular reflection of small objects 
on the ground. The pixels have proven to be useful in calcula- 
ting interband registration coefficients for all CCT formats. 

The coefficients indicate that the bands are registered to each 
other within a small fraction of a pixel. 

Author recommendations call for 1) better documentation, 
including that on tape files; 2) better calibration procedures 
and calibrated images; 3) nearest-neighbor resampling option 
in addition to cubic convolution; 4) computer processing using 
real-number formats tor all intermediate products, to eliminate 
the effects of integer roundoff; and 5) the use of video laser 
disk media for recording archival and browse images. Areas for 
future research are indicated. 


2 


KEYWORDS : 


Remote Sensing, Landsat-4, Thematic Mapper, Digital 
Image Processing, Calibration, Geometric Correction, 
Interband Registration, Computer Compatible Tapes, 
Data Formats, CCT-BT, CCT-AT, CCT-PT. 


T 


CONTENTS 




SECTION 


Page 


1 . 0 INTRODUCTION 

2.0 DIGITAL IMAGE PROCESSING f . . 

, 2.1 Raw Digital Image (CCT-BT) . . . . , . . 

2*1.1 Blocking Factor , 

2.1.2 Scan Shift . . * . 

2.1.3 Defective Detectors ...... 

2.1.4 Band-6 Format , 

2.1.5 Ancillary Data 

2.2 RadiomeLrically Calibrated Image (CCT-AT) 

2.2.1 Blocking Factor 

2.2.2 Scan Shift . 

2.2.3 Defective Detectors ...... 

2.2.4 Band-6 Format 

2.2.5, Calibration 

* 

2.2.6 Ancillary Data 

2.3 Geoco:;rected Image (CCT-PT) 

2.3.1 Blocking Factor 

2.3.2 Geometric Corrections Applied 

2.3.3 Earth Rotation Correction . . . 


9 

10 

12 

12 

13 

13 

14 
14 
16 
17 
17 
17 
17 


20 

20 

21 


3.0 RESULTS 

3.1 Specular Reflectance 

3.2 Registration 

3.2.1 Prime Focal Plane* Internal ...... 

3.2.2 Cold Focal Plane, Internal . . . . . 

3.2.3 Prime-Cold Comparisons 


22 

22 

22 

23 

23 

24 


4 


3.3 Striping 24 

3.4 CCT-AT Radiometric Abnormalities . , 25 

3.5 CCT-PT Radiometric Abnormalities , , 26 

3.6 Pseudocolor Pixel Prints 26 

4.0 DISCUSSION 27 

4.1 Summary , 27 

4.2 Research Areas ..... 27 

4.2.1 Overshoot and Undershoot 28 

4.2.2 Bright-Spot Recovery * . . 28 

4.2.3 Light Leaks 28 

4.2.4 Out-of-Band Spectral Response 29 

4.2.5 Specifications ............ 29 

4.2.6 Odd-Even Difference 29 

4.2.7 Within-Scan Droop 29 

4.2.8 Spectral Calibration 30 

4.2.9 Autocorrelation 30 

5.0 CONCLUSIONS 31 

5.1 Definitely True Conclusions 31 

5.2 Probably True Conclusions 31 

5.3 Possibly True Conclusions 32 

6.0 RECOMMENDATIONS 3 3 

6.1 Spacecraft Management 33 

6.2 Landsat-4 SCROUNGE 33 

6.3 Landsat-4 TIPS 33 

6.4 Landsat-4 NO A A 34 

6.5 Landsat-D’ 34 


7.0 REFERENCES CITED 


36 


ILLUSTRATIONS 


All figures follow the text, starting at page 40. 
FIGURES 


Figure 2-1. File organization of TM data and imager ^n the two 
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Figure 2-3. Ah example of data and organization of TM images in 
blocked CCT-BT format. 

Figure 2-4. An example of a Band 5 image in CCT-AT format. 
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Figure 2-5. An example of a Band 5 image in A* format. 

Figure 2-6. An example of a Pseudocolor Pixel Print image for TM 

data in B' format. 


Figure 2-7. An example of data for TM band 5, showing the 
defective channel #3. 

Figure 2-8. An example of data and organization for TM band 6, v 
in unblocked CCT-BT format. 


Figure 2-9. An Example of a Pseudocolor Pixel Print image for TM 
Band 6 in B' format. 


Figure 2-10. File organization of TM data and images on the two 
CCT-AT tapes used by LAS for data archive. 


Figure 2-11. Image file contents for the CCT-AT. 


Figure 2-12. An example of data and ogranization fur TM images 
in unblocked CCT-AT format. 

Figure 2-13. An example of TM data for spectral band 4 in CCT-AT 
format, unblocked, nominal calibration. 

Figure 2-14. An example of TM data for spectral band 2 in CCT-AT 
format, unblocked, nominal calibration. 


Figure 2-15. An example of TM data for spectral band 5 in CCT-AT 
format, unblocked, nominal calibration. 

Figure 2-16. An example of TM data for spectral band 6 in CCT-AT 
format, unblocked, nominal calibration. 
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Figure 2-17. Diagram of LAS calibration options for TM data. 

Figure 2-18. An example of TM data for spectral band 2 in CCT-AT 
tormat, unblocked, internal calibration. 

Figure 2-19. An example of TM data in CCT-AT format, unblocked, 
internal calibration. 


Figure 2-20. An example of TM data for spectral band 5 in CCT-AT 
format, unblocked, internal calibration. 
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Figure 2-21. An example of TM data for spectral band 6 in CCT-AT 
format, unblocked, internal calibration. 

Figure 2-22. The difference image between TM band 2 data 
processed using two different calibration options. 

Figure 2-23. The difference image between TM band 4 data 
processed using two different calibration options. 

Figure 2-24. The difference image between TM band 7 data 
processed using two different calibration options. 

Figure 2-25. The difference image between TM band 5 data 
processed using two different calibration options. 

Figure 2-26. Histograms showing the frequency of differences 
between ’mages produced by two different TM calibration options. 

Figure 2-27. Histograms for each channel showing the different 
rrequoncy of differences between band-7 images produced by two 
different calibration options. 


Figure 2-2,8. 
three CCT-PT 

Figure 2-29. 
format. 

Figure 3-1. 

Figure 3-2. 

Figure 3-3. 


File organization of TM data and images on the 
tapes produced by LAS as final products. 

Image file contents of Landsat-4 TM data in CCT-PT 

An example of spectral reflectivity. 

Interband registration for the primary focal plane. 
Interband registration for a 7-band TM image. 


Figure 3-4. An example of a histogram for one data channel 
showing empty bins. 


Figure 3-5. An example of empty bin locations and periodicity 
for data from detector 1, all 7 TM spectral bands. 
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SECTION 1 


INTRODUCTION 


The Landsat Thematic Mapper (TM) remote sensing team 
includes people with varied fields of interest and experience. 
There are the engineers, the scientists, the managers and the 
users, to name a few of the major groups. From a different 
perspective, there are the people who designed the system, those 
who built it (Covault, 1902a), those who tested it to see who the 
hardware and software lived up to expectations, and those who 
want to use the resulting images (Covault, 1902b; 1903) and 
digital products to solve problems (Williams, 1903; Everett et 
al., 1903; Abrams et al., 1983; and Quattrochi, 1903). The first 
three groups talk common languages in engineering, mathematics 
and computer programming. The fourth group, the users, in many 
cases may not be acquainted with or understand the technical 
documentation produced by the engineers, scientists and computer 
professionals . 

This document provides a transition. It is intended to 
provide a scientific characterization of the images, and infor- 
mation on the means by which the various images are being 
produced for the benefit of TM investigators, especially those 
who are concerned with radiometry. It assists users in infor- 
mation extraction by presenting some of the features and charac- 
teristics of the various digital image products used to produce 
the generally available final digitalimage product, the CCT-PT. 

It points to more detailed treatments of image characteristics or 
data form.ats in other documents. 

•k 

Evaluation of TM. image data will be the subject of another 
paper. Preliminary results by other investigators indicate that ' 
TM data have ”... an image resolution and quality beyond 
expectations, ... " (Covault, 1982b; 1983), due to performance at 
the optimum-expectation level of spacecraft and groundprocessing 
systems (S'alomonson & Koffler, 1983). Others have concluded 
"...the TM data exhibits very good data quality and reliability, 
and can be used for accurate mensuration studies ..." (Bernstein 
& Lotspiech, 1983) . 

Other papers in this symposium cover technical details of 
the Landsat-4 system (Barker, 1983). Calibration procedures are 
presented in Barker, Abrams, et al. (1984a). Absolute calibra- 
tion and engineering test results are in Barker, Abrams, et al. 
(1984b) while Barker, Ball, et al. (1984) describes post-launch 
calibration results for the TM reflective bands and studies 
improvements in the ground processing procedure. A discussion of 
sensor radiometry for the thermal band appears in Lansing & 

Barker (1984). Detailed discussion of results and recommen- 
dations is given in Barker (1984). This document will concen- 
trate on the characterization of digital products. 
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SECTION 2 — DIGITAL IMAGE PROCESSING 


Thematic Mapper (TM) in-orbit data is now available for 
general use (Covault, 1982b; 1983; Williams/ 1983), The SCROUNGE 
system, which comprised the Mission Management Facility (MMF), 
the Data Receive, Record, and Transmit Subsystem ( DRRTS ) , the 
Applications Development Data System (ADDS), and the Landsat-D 
Assessment System (LAS) has processed TM images on a research and 
development basis at Goddard Space Flight Center (GSFC), of the 
National Aeronautics and Space Administration (NASA) . SCROUNGE 
provided intermediate and/or final products directly to principal 
investigators (Pis) of TM assessment projects, and indirectly to 
the general public through the EROS Data Center (EDC) at Sioux 
Falls SD. SCROUNGE processing at NASA/GSFC provided images at 
the rate of one TM image per day. In July 1983, one year after 
launch, the TM Image Processing System (TIPS) took over image 
processing at a rate of 102 images per day. 

SCROUNGE-era computer-assisted (digital) image processing 
resulted in photoproducts and computer-compatible tapes (CCTs) 
representing different levels of processing. Most people will 
work with the images and CCTs at the final-product stage, but 
intermediate products are available as an engineering tool. TM 
data in final-product form has geometric and radiometric correc- 
tions applied and is ready to be used in image analysis; the 
intermediate products have different levels of corrections 
applied (Table 2-1). This section discusses the characteristics 
of the different types of digital products. 

In addition to the several high-density ( HDT ) tape formats 
used internally by NASA (Jai, 1982), SCROUNGE-era digital TM data 
exist in three formats (Table 2-1). The three formats are: 

o Raw rectified data (CCT-BT) ( Anon. , 1982a) is generally 
used only for internal transportation of digital data 
from one ground processing system to another, 
o Calibrated data ( CCT- AT) ( Baker & Gilson, 1982) should 
be used only by investigators doing radiometric 
characterization. 

o Geometrically corrected data (CCT-PT) (Baker & Gilson, 
1982) is the final product. The final product, the 
CCT-PT, is the one to which all of the various radio- 
metric and geometric corrections have been applied; 
this is the product that is ideal for image analysis 
for applications. 


10 


S>£ E 


» 
A** 

o 

3 

"W 

o 

hi 

J 

C! 

F-i 

<*3 

IKK 

« Q 

« * 
o 
& 
a 

C3 

«-i 

o o 

y lPiLI 

— g 
Q £ 
t* 5 
(?) g 
K »s 
(?) I* 

(3 

w 

"0 

c 

<0 

•J 


© 

w 

■ EH 

2Si 

w 

©. 


CO 

CO 

05 


Csl . 

s£ 
-< 
W Q 

mJ 

< < 

e5 tl 

§2 

LU O 

ce h 

LU ' 

I “ 

UJ d 

gS 

Sr W 

3 Q 

oS 

w 

h- 

O 

o 


0 . 

►" 

o 

u 


h* 

03 

(j 

O 


(O 


a 

n 

o 

in 

CM 

to 


a 

o> 

CD 


M 

O 


(S 

ca 


CM 


a 

A 

o 

Ln 

CM 

<C 


CO 

Ui 

CL 

IS 


03 

2 

lit 

a 


H 

< 

£ 

cc 

o 

LL 

a 

2 

« 

ca 


ui 

•ml 

LL 

*•'. 

03 

Q 

a 

o 

a 

UI 

cc 

111 

o 

< 

£ 


03 

UJ 

s 

ca 

R 

CO 

CO 

CM 






03 

h 




UI 

L. 

< 

IM 



fc 

£ 

a 



ca 

o 

u 

A 

S? 

a 

03 

ca 

rs 

05 

s 


CM CD 

V 

r- 

CM 


03 


co 

CD 


C 3 

05 

CD 


CD 

in 

CO 

CO 

£ 

5 


Q 

in 


X 

h- 

O 

2 

ui 

mJ 

a 

cc 

o 

o 

UJ 

cc 


X 
H 

a 
2 2 
UJ c 
aJ 


CC 

g 

o 

2 

a 


a a 

UJ UJ 
-1 mJ 

cl a. 
£ £ 
< < 
03 03 
UJ UI 

cc cc 


UJ 

03 g 
uj O 
> 2 


03 

2 

< 

a 

03 

UJ 

03 


UI 

2 


a 

o 

-j 

ca 


a 

2 

03 

03 

UI 

O 

o 

cc 

CL 

UI 

CC 

CL 


CC 

tu 

2 

o 

5 

< 

CC 


o 

5 

tn 

s 

o 

ut 


03 

UJ 

> 


03 

UI 

> 



1 

o 


o 

o 


UI 

> 



o 

£2 

ca 

s 

UJ 

UJ 

CC 


d 

£ 

u 

o 

o 

2 

CL 


03 

ca 

05 

T— 

r- 

o 

CM 

CO 

o 

S in 

O 

2 

3 

u. 

o 

2 


CC 
UJ 
H 
2 
uj 

<5 

»- d 
<C ca 
a < 
o)=i 

« g| 

uj < 


M 

o 

£1 

rtt 


General Comparison among Landsat~4 Thematic Mapper Digital Products. 


2,1 Raw Digital Image (CCT-BT) 

# , _ , , . „ 


SCROUNGE-era partially processed TM data was produced in a 
format (CCT-BT) used for data exchange between the ADDS and LAS 
(Anon./ 1982a). CCT-DT tapes are not ready for general applica- 
tions digital image processing. They are available through the 
NASA Science Office only to TM principal investigators (Pis), 
Radiometric correction and geometric registration procedures have 
not been applied to the image data. Start-line pulses and 
initial image pixels, not present in all channels of the tele- 
metry data, have been trimmed before writing the image t,o the 
CCT-BT, Reverse scans have been flipped so that east and west 
ends of each scan are properly oriented. The interleaved 
minor-frame format used in the telemetry was sorted into 
band-sequential format (Figure 2-1), 

The images were written as 8-bit byte data on two standard 
6250 bpi CCTs in band sequential (BSQ) format. Four image files 
are on the first CCT-BT, and three image files are on the second 
(Figure 2-1). Tape identification files are also placed on both 
tapes. In addition, the first tape contains a special data file 
with housekeeping, ephemeris, attitude, and minor scan data plus 
data quality information (Section 2.1.5). 

Each image line follows a particular format (Anor.., 1982a). 
The first 10 bytes of each imago line provide support data, 
including the scan number for that TM scene. The next 6176 
pixels are image data, with the first pixel at the western edge 
of the scene. The image data are followed by 64 bytes of zero 
fill, by 64 bytes of support data containing channel gain and 
off set . information , and by a final 64 bytes of zero fill (Figure 
2-2). This yields an image line of 6400 bytes. Five such lines 
make a record on the tape (Section 2.1.1). 


2.1.1 Blocking Factor 

To save space on the tapes, image lines were not written as 
separate records to tape files, but were blocked five (5) lines 
per tape record. If an image is copied fro;,* the tape without 
taking into consideration the blocking factor, the image is 
compressed by having only one line out of five in the vertical 
direction and five lines across in ( the horizontal direction 
(Figure 2-3). The image has a line length (logical record length 
= LRECL in IBM terminology) of 6400 bytes; the physica- record • 
length (« BLKSIZE ) on the tape is 32000 bytes. Some mini- 
computers are u'nable. to read blocks of this size. 

The image must be unblocked to be useful (Figure 2-2). 

The unblocking may be a preprocessing step on an auxiliary 
computer where the image is copied from one tape to another by a 
utility program or unblocking may be performed automatically by 
the computer program that puts the image or a portion of it on a 
disk file for computer access. 
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2.1,2 Scan Shift 

TM imago data in B format is not aligned between forward 
scans and reverse scans (Engel, 1980). A linear feature is 
> -operly recorded within each scan, but between scons there is a 
large offset (Kieffer, Abrams, et al,, 1983) (Figure 2-4). 

The offset of the image in B format can be corrected to 
yield a nominal alignment of pixel columns, a B’ format (Figures 
2-5 and 2-6). To do this, an analyst makes repeated runs of vt 
"mosaic" program; one run for each 16 image linfcs. The number of 
pixels offset is variable, dependent upon scene latitude. An 
offset that appears ideal for Mississippi does not appear to be 
appropriate for Maryland nor for New Mexico, 

•The offset applied to the B format image corrects the 
alignment between forescan and backsean, but does not take into 
account all factors. It is necessary to correct the B format 
image for the > station of the earth. The correction for the 
rotation of the earth is applied by moving the start sample for 
each mosaic of 16 lines to the east (i t e., by increasing the 
value of the start sample by a constant). The value of the 
constant used is a function of angular distance from the equator. 
Operationally, it is easiest to produce an intermediate image, to 
determine the correction for the rotation of the earth, and then 
to apply that correction to the image to be studied. The appli- 
cation of both scan-offset and earth-rotation corrections 
produces a B M format image. 

The image in B format is available through the NASA Science 
Office to TM Pis only. The image in B' and B* ' formats must be 
picvd.ited by the individual investigator doing his own digital 
image processing? the B' and B 1 ' image formats are not standard 
Lundsat-4 digital products. 

Note that the B' ' image format is close to being 
geometrically correct, without any resampling. It is not a 
precision-registered, geometrically-correct image. This format 
is useful for radiometric studies that require ground truth. It 
is also useful for other studies where lack of ancillary data 
prohibit the construction of a final-product image. 


2.1.3 Defective Channels 

Before launch, four channels were found to be defective. 
Because of construction difficulties, the sensor was launched 
with the known defects, and ground-processing software was 
adjusted to replace defective data. In the B format image, these 
corrections have not been applied. 


Detector 3 of band 5 does not respond to radiance input 
(Barker, Ball, et al., 1984). Instead of valid data, the 
instrument counts or data numbers (DNs) in that data channel are 
at least one order of magnitude lower than that recorded in the 
previous channel and the following channel (Figure 2-7). ‘ 
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Channel 2 of band 2 and channel 7 of band 7 have twice as 
much noise as the average for all channels for the respective 
band (Barker, 1983 and Barker, Ball, et al,, 1984), Channel 4 of 
band 2 has a smaller modular transfer function (MTF) than desired 
(Barker, 1983 and Barker, Ball, et al., 1904), 


2,1,4 Band~6 Format 

Band 6 is the thermal band, with a larger instantaneous 
field of view (IFOV) than the ref lected-light bands (Engel, 1900; 
Covault, 1982a), The "footprint" for band 6 is 120m by 120m 
rather than 30m by 30m; this results in an array of numbers for 
the image that, is 1/16 the size of the array for the other bands. 
For digital image processing purposes, the array size of each 
band should be ,the same, To accomplish this, processing software 
replicates each band-6 pixel to produce four pixels along a line 
and four identical lines. 

In the B format image, the replication process is only half 
done. The pixel values are replicated in the scan direction, but 
the scan line has not been replicated (Figure 2-8). As a result, 
in the B format data, the band 6 image file has only 1/4 as many 
records as the other files (Figure 2-1). In the B' (Figure 2-9) 
image format, the scan ‘lines have been duplicated by the user so 
that the band-6 file has the same number of records. 


2.1.5 Ancillary Data 

* 

In addition to the image files that are characterized in the 
above sections, each tape produced in SCROUNGE-era B format 
contains data files that have useful ancillary data. . 

The first ancillary data file is the TAPE ID file (Anon., 
1982a), found as the first file on both tapes. The ID file is 
one physical record, or five logical records (Section 2.1.1). A 
detailed description of the information and format of the ID file 
appears in Table 2-2. 

The second ancillary data file is the Header Ancillary 
Annotation Trailer (HAAT) file, found as the second file on the 
first tape only (Anon., 1982a). The HAAT file is seven physical 
blocks in length? each logical record is 6400 bytes. HAAT file 
internal head and tail records contain general timing and loca- 
tion information; other records contain detailed housekeeping, 
ephemeris and attitude data. The general information is 
described below; characterization of the detailed information can 
be found in the "Interface Control Document" (Anon., 1982a). 
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Table 2-2. File Contents for the TAPE ID file. 


Label 

Location 

Example 

Note 


(bytes) 



SYSTEM ID & DATE 

1-12 

L4TCC8230801 

1 

IMAGE FORMAT & VOLUME 

15-22 

CCT-BT 1 

2 

IMAGE SCENE ID. 

23-38 

0109 015-033 

3 

SOFTWARE VERSION 

41-52 

ADDS0 102.00 



Notes: 1. Landsat-4, Thematic Mapper, Computer Compatible Tape 

created on day 308 (4 November) of 1982, tape #1. 

2. First tape of data in B format, 

3. Image taken on Launch Day + 109 days, for World 


Reference System (WRS) Path=015, Row=033. 


The first-, record in the HAAT file contains general timing 
information as literal or alphanumeric fields. The times at 
which the first scan and the last scan lines for this scene were 
taken are given in bytes 21-36 and 37-52 in YYDDDHIIMMSSTTTFF 
format. In this format "YY" is the year, ignoring the century 
numbers, "DDD" is the Julian day of the year, "HH" is the hour of 
the day and "MM" is the minute of the hour. The seconds are 
reported by "SS" and the thousandths of a second by TTT " ; 

"FF" values are the 16ths of a millisecond (Jai, 1982). The 
start and stop Payload Correction Data (PCD) t.ijnes are' given in 
bytes 53-68 and 69-84 in the same format. 

The seventh or last record in the HAAT file contains general 
location information as literal fields. The date of acquisition 
is given in bytes 5-11. The latitude and longitude for the image 
center are given in bytes 13-28, and the same for the MW corner 
in bytes 39-54. The World Reference System (WRS) path-row coor- 
dinates are given in bytes 31-37. Bytes 66-78 contain solar 
elevation data. The full scene identification code may be found 
in bytes 92-117. 


2.2 Rad iometrically Calibrated Image (CCT-AT) 


The intermediate digital-data product during the 
SCROUMGE-era was the CCT-AT, In producing the CCT-AT digital 
product, some of the imperfections of the CCT-BT images are 
removed, and the images are written on the tapes in a general 
user format (Anon., 1979). The images are not yet suitable for 
general applications analyses. 

Radiometric correction procedures have been applied to 
the image data (Barker, Abrams, et al . , 1984a), but not geometric 
registration procedures. CCT-AT tapes, which are generated by 
NASA-GSFC for archival purposes, are available through EDC to 
researchers interested in radiometry. 

The images were written in band sequential format ( BSQ) 
as 8-bit byte data on two standard 6250 bpi CCTs. Three image 
files are on the first CCT-AT, and four image files are on the 
second (Figure 2-10). Tape identification files and other data 
files are also placed on both tapes. Identification files, data 
files and image files are packaged in the manner recommended by 
the Landsat Ground Station Operators Working Group (LGSOWG) 
(Anon., 1979; Baker, et al . , 1983). 

The formats for the image file and for the image lines 
are different for a CCT-BT (Section 2.1) and a CCT-AT (Figure 
2-11). The image file for a CCT-AT has prefixed one internal 
header record to the image data. In the A format image, the 
image lines do not have the initial 10 bytes of labeling infor- 
mation, but they do have an appendage of 480 bytes (Section 
2.2.1); the first image pixel is at byte #1. 


2.2.1 Blocking Factor 


CCT-AT image data during the SCROUNGE-era were blocked 4 
lines per tape record {Baker, et al., 1903). In addition, 
hardware efficiency demands in LAS require that the 6176 data 
pixels per line be padded so that the line length is an even 
multiple of 512. This results in an additional 480 bytes per 
line, for a total line length of 6656 bytes and a physical record 
length (= BLKSIZE) of 26624 bytes (Figure 2-11). 

The one-line internal header record for the image file 
is in a separate record the same length as the image line, Most 
of that record is blank. When the internal header record is 
reblocked along with the true image lines, it produces one line 
with alphanumeric header data and three lines with "spaces" 

(ASCII code = 32) (Figure 2-12). 

The- resulting image array does not have its start line and 
start sample at location 1,1 but rather at line=5, sample^l. 


2.2.2 Scan Shift 

The scan shift is the same in both A format and B format 
images. Since the transformation from B to A does not involve 
geometric corrections (Figures 2-13 to 2-16), A' and A' 1 
products are derived by the user from an A format image in the 
same way as B‘ and B’ 1 images were derived from a B format image 
( Section 2.1.2). 


2,2.3 Defective Detectors 

In the A format image, a software "quick fix" was applied 
by SCROUNGE to two of the channels of data produced by defective 
detectors (Section 2.1.3). Data from band 5, channel 3, is 
replaced with data from channel 4, and data from band 2, channel 
4, is replaced with data from channel 5. This results in two 
identical lines out of every 16, in both forward and reverse 
scans in images from bands 2 and 5 (Figures 2-14 and 2-15). 


2.2.4 Band-6 Format' 

Unlike the band 6 image in B format, the band 6 image file 
in A format has the same structure as the other bands (Figures 
2-10 and 2-16). The size of the image array was enlarged to 
provide array compatibility with that of the other bands. 


2.2.5 Calibration 

TM image data in A format is produced by radiometric 
calibration (Barker, Abrams, et al., 1984a) of images in B 
format. 


Prior to launch, the channels were calibrated against a 
National Bureau of Standards integrating sphere {Engel, 1980; 
Anon., 1982b; Barker, Abrams, et al., 1984a; Barker, Ball, et 
al., 1984). The calibrated channels were used to calibrate the 
on-board internal Calibration (IC) system (Barker, Abrams, et 
al., 1984a; Barker, Ball, et al., 1984). Since launch, data from 
the IC was used to determine gain and offset for each 
channel (Barker, Abrams, et al., 1984a; 1984b). The application 
of sensor-specific calibration factors provided by the IC to the 
reflectance data allows production of high-quality images 
containing accurate radiometric data (Barker, Abrams, et al . , 
1984a). The actual application of calibration data proceeds at 
high speed through the use of a radiometric look up table (RLUT) . 

A B-format image can be calibrated in several ways 
(Figure 2-17) to produce an A-format image. For example, the 
authors placed a special order for a particular image to be cali- 
brated by using pre-launch nominal values, without use of the IC, 
plus histogram equalization procedures to produce one set 
(Figures 2-1,3 to 2-16) of A format images, and by the IC plus 
histogram equalization procedures to produce another set (Figures 
2-18 to 2-21) of A format images. An investigator interested in 
radiometry can make special requests for image calibration in A 
format : 

o An image that has received nominal calibration only; 

this requires special processing by the ADDS portion of 
the SCROUNGE system. 

o An image that has received IC calibration only, without 
histogram normalization. 

o An image that has received the nominal calibration, 
followed by histogram equalization; this requires 
special processing by the ADDS portion of the SCROUNGE 
system. 

o An image that has received IC calibration, followed by 
histogram equalization. Under the SCROUNGE system, 
this is the option used if there is no special request. 

o An image that has received histogram equalization 
processing only. 

The authors found differences between the products of the 
two types of calibration, IC followed by histogram equalization 
and nominal followed by histogram equalization. These 
differences were investigated by producing difference images for 
corresponding bands between the IC calibration and nominal cali- 
bration A format images. Examples -of the difference images 
(Figures 2-22 to 2-25) show that: 1) differences between the two 

calibration techniques are almost always less than 2 DN levels 
(see below); 2) changes in calibration tend to be autocorrelated 
along scan lines, since the corrections have the same calibration 
offset for all data produced by a given channel, and since the 
calibration gain for similar pixels is similar; 3) each scan line 
tends to have changes in calibration that are independent of what 
happens in neighboring lines; and 4) each TM band has its own 
pattern of changes in calibration. 
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Histograms have been produced for calibration difference 
images for all TM bands (Figure 2-26). These histograms show 
that for all TM bands except band 4, most pixel DNs are not 
changed by the change in calibration procedure. TM band 4 has a 
1-DN change between the two calibration procedures for 86.5% of 
the sampled image pixels and a 2-DN change of 5.2%, This could 
arise from the unusual IC pulse shape for band 4 channels (Barker 
et al, 1983b). TM band 1 has a 1-DN change for 33.2% of the 
pixels. TM band 3 has a 1-DN change between the two calibration 
procedures for 16,4% of the image pixels. Band 2 has a 1-DN 
change for 15.4% and band 5 has a 1-DN change for 8.7% of the 
pixels. Band 7 has 2-DN changes for 0.1%. Other DN changes 
involve trivial percentages. 

Histograms were also produced for each of the 16 channels 
for each spectral band (Figure 2-27). These show that the two 
types of calibration produce different results for each channel. 

A word of caution. Instrument variability with time 
(Barker, Abrams, et al,, 1984b) has placed constraints on the 
several calibration procedure options available (Figure 2-17). 
Radiometric image analysis (Barker, 1983; Barker, Abrams, et al., 
1984b) has indicated that it is essential to apply IC radiometric 
correction before other calibration corrections to preserve the 
radiometric accuracy needed for applications analyses. 


2.2.6 Ancillary Data 

In addition to the image files that are characterized in the 
above sections, each tape in A format contains v data files that 
are not, in image format. These data files contain information 
that may be necessary or very useful to the individual interested 
in information extraction. 

The first ancillary data file is the LGSOWG VOLUME DIRECTORY 
file (Anon,, 1979; Baker & Gilson, 1982), found as the first file 
on both tapes. The file has 17 physical records, or 68 logical 
records (Section 2.2.1). The first logical record contains the 
LGSWOG Volume Descriptor Record. The next 16 records contain 
pointers to various files on that tape. Details of file contents 
may be obtained from Baker & Gilson (1982). 

The second ancillary file is the LAS label file for HAAT 
data. This label file contains a record describing the HAAT data 
file and optional history records. 

The third ancillary data file is the HAAT file, found 
as the third file on the first tape only (Baker & Gilson, 1982). 
HAAT file records contain detailed housekeeping, ephemeris and 
attitude data. The general information is found on records #2 
and #32, which are copies of records #1 and #7 from the HAAT file 
of the CCT-BT tape (Section 2.1.5). Characterization of the 
detailed information may be obtained in Baker & Gilson (1982). 
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The fourth ancillary data file is the LAS label file for 
image data, band 1. The general contents of the label file are 
the same as that described above for the HAAT label file. 

Files with this same format and data are found as external 
labels in front of each image file on both tapes (Figure 2-10). 


2.3 Geocorrected Image (CCT-PT) 


The final and finished digital-data product is the CCT-PT. 
The radiometrically fully-corrected images were written on the 
tapes in a user format (Anon., 1979) and are suitable for general 
applications analyses. 

Completely processed TM data generated by SCROUNGE (Baker & 
Gilson, 1982) is available in the CCT-PT format used for product 
data. Radiometric and geometric correction procedures have been 
applied to the image data. This is the format for the image data 
on tape that NASA/GSFC supplied during the SCROUNGE era to the 
EROS Data Center ( EDC ) . 

The images were written as 8-bit byte data in BSQ format on 
three standard 6250 bpi CCTs (Figure 2-28). Two image files 
are on the first and second CCT-PT, and three image files are on 
the third. Tape identification files and other data files are 
also placed on the tapes. Identification files, data files and 
image files are packaged in the manner recommended by LGSOWG (see 
Baker & Gilson, 19S2a) . 


2.3.1 Blocking Factor 

To save space on the tapes, image lines were not written 
separately to tape files, but were blocked 4 lines per tape 
record (Baker & Gilson, 1982). In addition, hardware efficiency 
demands within the LAS computer required that the 6967 data 
pixels per line be padded so that the line length was an even 
multiple of 512. This padding resulted in an additional 201 
bytes per line, for a total line length of 7168 bytes and a 
blocksize of 28672 bytes (Figure 2-29). 

The one-line internal header record for the image file is in 
a separate record in front of the image data. The situation 
is nearly identical to that for A format images (Section 2.2.1). 

The resulting image array does not have its start line and 
start sample at location 1,1 but rather at line=5, byte=l. 


2.3.2 Geometric Corrections Applied 

Geometric correction places TM image samples into a map 
coordinate system. This simplifies subsequent applications data 
processing (Beyer, MS), but does not affect the spectral response 
of ground cover types (Bartolucci et al., 1983). 
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During the SCROUNGE-era , ground processing produced 
geometrically corrected images using resampling based on syste- 
matic correction data (SCD), derived from payload correction data 
(PCD) (Anon., 1902b) and mirror scan correction data (MSCD) . 
During this processing, flight segment telemetry was converted 
into a geometric projection of the scan pattern (Beyer, 1983). A 
transformation was generated between the geoid location of each 
pixel and the appropriate coordinates for any output grid system. 
The output grid was filled with data using a resampling tech- 
nique. SCROUNGE used cubic convolution procedures rather than 
nearest neighbor or bilinear interpolation for resampling. The 
original 30 m IFOV of the pixels in the B and A format images was 
resampled to a grid interval of 28.5 m. 

Precision geometric correction using the known locations of 
ground control points was not available during the SCROUNGE-era . 
The nominal corrections produced from satellite ephemeris data 
provide images that closely correspond to maps. Working with a 
CCT-PT image of Death Valley CA, registered to a topographic 
base, Abrams et al. (1983) found a registration error of less 
than two pixels, for terraine with 7000 feet of local relief. 
McDonald (1983), working with a CCT-PT image of flat farmland in 
Arkansas, reported an error in geometric accuracy of about 0.1 
pixel) Goodenough (1983) found that ground registration errors 
along a scan were greater than across a scan (32m versus 40m). 
Bryant et al . ( 1983) calculated an average vector offset of 
132.77m using 75 ground control points for a full scene. 

Some apparent errors in geometric precision in film products 
result from imperfect geometric fidelity of filmwriters. Kiefer 
et al. (1983) found scale errors corresponding to 4 pixels and 
skew errors of 3 pixels. Large errors reported (68~72m) were 
attributed at least in part to "highly significant skewness" in 
the rotary-drum, film-writing device (Bender, et al . , 1983). 


2.3.3 Earth Rotation Correction 

• A consequence of applying geometric corrections is the 
correction for the rotation of the earth under the orbiting 
satellite. This correction results in the two black wedges found 
in full-scene or near-edge photoproducts and the corresponding 
fields of zero values found in the digital image files as "fill" 
pixels (Figure 2-29). These two fill fields have caused problems 
in image enhancement procedures that expand the dynamic range of 
raw or transformed data by thresholding selected percentages at 
the high and low ends of the frequency distribution. 

Because the effects of the earth's rotation are corrected 
by using fill pixels, the first image pixel will be found in tape 
record line 5, but at a sample number some distance in from the 
start of the line. The distance is a variable depending upon 
latitude north or south. If it is necessary to locate the first 
image pixel on any line, the analyst or a computer program must 
search for the first non-zero pixel in that line, and then move 
three pixels to the right, to avoid edge effects caused by cubic 
convolution resampling of image data and fill data (Baker & 
Gilson, 1982). 
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SECTION 2 — RESULTS 


The authors employed a number of digital image analysis 
tools to examine particular features of Landsat-4 TM images. 
Preliminary results from these investigations are presented 
below. 


3.1 Specular Reflectance 


Statistically-based image evaluation procedures commonly 
show the presence of a few very high values, often reaching 
saturation levels (DN=255). Line and sample location data and 
printouts of pixel values show (Figure 3-1) that the same pixels 
are involved in each band. This indicates that the high DN 
values result from specular reflectance from some object within 
the scene and are not noise. 

In a few cases, the source of specular reflectance was 
determined. In a Goddard Space Flight Center (GSFC) TM scene, 
curved surfaces (the inflated warehouse near Bldg. 20, the water 
tower at Bldg. 16, and the pressure dome for Bldg. 15) show 
specular reflectance. In the same scene, some flat roof areas 
also show specular reflectance. Specular reflectance could be 
obtained from silo tops, from covered stadiums, from cars parked 
v in alignment in a parking lot, from pitched roofs, etc. 

The specular-reflectance pixele provided an opportunity to 
look for instrument lag. It is important to know if the TM 
sensors are slow to recover from a bright reflectance reading. 
Examination of the DN values around several specular-reflectance 
pixels indicates that the next pixel in the scan line is at least 
100 DN values lower. Any adjacent pixel with DN values above 
those of its neighboring pixels can easily be explained as a 
mixed pixel from a large bright source. Statistical studies of 
large numbers of specular-reflectance pixels are being 
considered, but at present there is no clear evidence for 
instrument lag. 


3.2 Registration 


Two aspects of registration can be considered: 1) the 

band-to-band registration that provides an accurate multispectral 
image array; and 2) the image to ground registration that allows 
the image to be used as a map. The latter aspect was not inves- 
tigated by the authors; it is discussed in Section 2.3.2. 
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Investigation o f specular radiometric responses (Section 
3.1) has shown that there is close interband registration 
in both the B and A format images, which have the same image 
geometry, and in the P format images. Identification of specular 
reflectance pixels is complicated by the different sensitivities 
of the different spectral bands (see also Bryant et al., 1903a), 
In some bands, notably bands 1 and 4, the spectral reflectance 
pixels are saturated; saturated pixels record a maximum DN of 
255. In other bands, notably 2 and 3, the spectral reflectance 
provides a peak DN value considerably above that of its neigh- 
bors, but not at saturation. The same peak can almost always be 
identified in all spectral bands (Figures 3-2 and 3-3). 


3.2.1 Prime Focal Plane, Internal 

The two early, 4-band TM scenes for Detroit, MI (scene 
ID * 40004-15401 and 40009-15413) in CCT-BT format, and one 
7-band scene for Memphis, TN (ID = 40037-16033) in CCT-AT format 
have have been examined for interband registration as indicated 
by spectral reflectance. Each scene was treated separately. 

In some cases, specular reflectance pixels for large areas such 
as quarries occupy too large and diverse an area for registration 
evaluation. In most cases, specular reflectance involves only 
one or two pixels, and these locations are suitable for interband 
registration evaluation. 

The authors have made center-of-gravity calculations for 
small image segments centered on isolated specular reflectance 
pixels. The assumption is that the weighted-average line and 
sample values for each band in the multispectral small image 
segment will have the same value if the bands are corrected. 

The data for each Detroit scene taken separately indicates that 
the bands on the primary focal plane (Engel, 1980) are registered 
to a fraction of a pixel in both line and sample directions 
(Figure 3-2). The Detroit-1 image had a registration error of 
up to 0.14, pixels in the scan direction and up to 0.28 pixels 
across scan; the equivalent values for the Detroit-2 image are 
0.27 and 0.48. Errors of up to 0.45 pixels in scan and up to 
0.15 pixels across scan were obtained for a comparison of bands 
1-4 for the Memphis scene (Figure 3-3). 

Other authors have reported similar results. Using a 
different procedure, that of line-to-line displacement, Fusco and 
Mehl (1983) found a misregistration that measured less than 0.1 
pixel. The same value was determined by Bryant et al. (1983a; 
1983b). Visual comparison of contrasting targets using an inter- 
active video display on a part of the Washington D.C. scene (ID = 
40109-15140) found that bands 1-4 were perfectly coregistered 
(Bender et al., 1983). 


3.2.2 Cold Focal Plane, Internal 

Center-of-gravity calculations for specular reflectance 
pixels for a subimage from the Memphis scene (ID = 40037-16033) 
in CCT-AT format indicates that bands 5 and 7 on the cold focal 
plane (Engel, 1980) are registered to 0.04 pixels in scan and 
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0,05 pixels across scan directions (Figure 3-3). Band 6 was not 
investigated at this time because it senses emitted rather than 
reflected light. 

Band-to-band registration error of bands 5 and 7 was 
reported to be less than 0.1 pixel (Bryant et al., 1983a; 1903b), 


3.2.3 Prime-Cold Comparisons 

Most studies indicate that the reflected-light bands on the 
two focal planes are only slightly misregistered . The authors' 
center-of-gravity (Section 3,2.1) comparison between the two 
focal planes using a part of the Memphis scene (.’)*= 40037- 
16033) in CCT-AT format indicates that the error -n registration 
of the reflected-light bands to each other is 0 . 1 pixels in the 
scan direction and 0,3 pixels across scan (Figure 3-3). Card et 
al. (1983), using a moving-block correlation method on the 
Arkansas scene (ID « 40037-16031), reported errors of 0.25 and 
0.50 pixels along and across scan. Using the phase correlation 
image alignment method on the 2 November 1982 Washington D.C. 
scene (ID = 40109-15140), Bryant et al. (1983a; 1983b) found 
between 0,1 and 0.35 pixels error in registration. 

Other authors, using different images and techniques, have 
reported larger registration errors. Fusco and Mehl (1983) 
reported a displacement of 0.5 pixels between focal planes for 
the Arkansas scene (ID = 40037-16031) , as measured by correlation 
profile maxima. Bender et al. (1983) found a displacement of 1-2 
pixels by interactive video display analysis of the Washington 
D. C .• scene (ID = 40109-15140). Barker ( 1983) repo.vted an 
apparent misregistration of about 0.75 pixels along scan and 0.2 
pixels across scan. Using an analysis of linear features, Gurney 
& Eng (1983) reported error values of 0.75 pixels and 0.5 pixels 
along and across scans respectively. 


3.3 Striping 


Digital image processing for image segments that are very 
uniform shows the presence of light and dark stripes (Metzler & 
Malila, 1983). Some images have forward and reverse scans that 
are consistantly different by up to 4 DNs (Barker, 1983; Metzler 
& Malila, 1983). Other images show only minor differences 
between forward and reverse scans; Kiefer et al. (1983) reported 
a difference of about 0.1 DN in the mean for a CCT-AT image 
segment from the Washington D.C. scene (ID = 40109-15140), and 
Bender et al. (1983) reported DN differences with a range of 
0.3-0. 9 for a different part of the same scene. 

Another type of striping is found within each scan. Some 
image lines have DNs that are consistantly higher or lower than 
those of neighboring lines for the same scan; Kiefer, Abrams, et 
al . (1983) found a typical 0.8 DN difference in the mean among 
the 16 channels in each band. 


In phofcopr Jducts, the range in DN values results in 
line-parallel stripes by line and by scan/ seen in water bodies/ 
ice and snow fields/ and other large uniform areas. Striping is 
present in 3-band color composite images/ 1-band pseudocolor 
images/ and black and white images. 

The sensor pecularities that produce the striping are also 
present in data from heterogenous areas/ but in these cases/ 
the image signal variability overpowers the striping noise. 

* 

Striping is evident in B, A, and P format images, Because 
of the difficulties of making sensors exactly alike, it is not 
suprising that striping is found in B format images. Striping in 
A format images indicates that sensor calibration is not yet able 
to remove all of the differences in sensitivity of neighboring 
detectors. Striping in P format images indicates that cubic 
convolution resampling does not remove the effect. 

Photointerpretation of the image data is only slightly 
hampered by the striping, but striping poses major problems for 
digital image processing. Image enhancement procedures will 
enhance the striping. Image classification mapping procedures 
will produce incorrect classifications/ most classification 
procedures work one pixel at a time and are unable to compensate 
for the radiometric changes that produce light or dark stripes. 


3.4 CCT-AT Radiometric Abnormalities 


One consequence of applying radiometric corrections to byte 
data is the phenomenon of empty bins in frequency histograms 
(Bartolucci et al . , 1983/ Bender et al,, 1983). A histogram of a 
large area in a TM scene (Figure 3-4) in A format has occasional 
DN values in the middle of a gradational sequence with no 
frequency counts. This is especially true if the histogram is 
otained for a single channel of data. 

This is an effect of integer arithmetic. Calibration 
corrections are not integer numbers, but numbers with decimal 
fractions. Image DN numbers are integers, ranging from zero to 
to 255. When an integer is multiplied by a decimal value, and 
the result is stored as an integer value, there will be certain 
DN values in the resulting calibrated image that can not be 
generated. A histogram of the calibrated data will have a zero 
tally for these particular DNs. 

The cubic convolution resampling procedure used to produce 
the final product image produces new DN values (Section 2.3.2), 
so that empty bins do not exist in P format images (Metzler & 
Malila, 1983). 

In a histogram for any data channel, the empty bins occur 
with a constant periodicity (Figure 3-5). The period of 
occurrence is different for each channel and each spectral band. 
The difference in the reoccurrence period arises from the 
different calibration factors applied to data from each detector 
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of each band. The periodicity resuits from the linear calibra- 
tion correction procedures used, 

The periodicity in the data for a given channel can be used 
to back-calculate one of the calibration factors used for that 
channel. For example, data channel 1 of band 7 for Memphis 
(scene ID ■ 40037-16033) has a period of 33 for the presence of 
empty bins. This means that a range of 32 DNs in the 
uncalibrated image is placed into a range of 33 DNs in the 
calibrated image. Thus 

33/32 * 1.03125 

is calculated to be the relative calibration gain factor, for 
original digital counts to calibrated digital counts. Neither 
the calibration offset factor nor the original radiance values 
can be determined in this manner. 


3,5 CCT-PT Radiometric Abnormalities 


An artifact is created by applying cubic convolution 
resampling methods to digital TM data. At each end of a scan of 
16 lines, the CCT-PT image has fill pixels with zero DNs (Suction 
2,3.3). The resampling procedure operates in such a manner that 
the near-edge combination of zero DNs and real-data DNs produces 
intermediate DNs along a 3-pixel border zone (Baker & Gilson, 
1982). These "fuzzy frame edges"- (Kiefer et al., 1983) result in 
erroneous histograms and misclassif ied pixels when the image is 
processed for applications analyses. 


3.6 Pseudocolor Pixel Prints 


The authors have devised a new tool for the radiometric 
evaluation of TM images that takes advantage of color graphics 
and human response to changes in color. The pseudocolor pixel 
print (PsCPP) product is produced by applying a custom color 
look up table to an image. Standard pseudocolor (Ps) look up 
tables span the range of 0-255 in a gradient that allows gross 
structure to be easily detected but that masks fine structure. 

The PsCPP look up table was designed to show differences in DN 
values of one count as color differences, over a repeating cycle 
of 32 colors. PsCPP uses eight cycles of colors to span the 
entire range of 0-255; this masks coarse structure, but the com- 
bination of Ps and PsCPP processing permits analysis of both 
levels of structure. 

PsCPP products show dramatically the fine structure in TM 
single-band images (Figures 2-6 and 2-9). Banding becomes very 
obvious, especially that caused by dead detectors (Section 2.1.3) 
rather than by minor variation in detector response and cali- 
bration (Section 3.2). 
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SECTION 4 -« DISCUSSION 


4.1 Summary 


Thematic mapper (TM) digital data has characteristics of 
tape files, image format and image radiometry that are Important 
to the applications scientist. First of all, the data resides in. 
different tape files at initial (Figure 2-1), intermediate 
(Figure 2-10), and final (Figure 2-28) stages of SCROUNGE 
processing at LAS. Within each image data file, the image lines 
are blocked by a factor of either 5 for a CCT-BT or 4 for a 
CCT-AT and CCT-PT (Sections 2.1.1, 2.2.1, and 2.3.1); access to a 
mainframe or large minicomputer is needed because of the large 
block sizes involved (i.e., up to 32,000 bytes). On each set of 
tapes, the image data files are in correct spectral sequence, but 
not in correct number sequence; band 7 is located before band 6 
because band 7 was added to the sensor package late in the 
program, and is thus the 7th set of sensors (Covault, 1982a) . 
Within each CCT format, the image file has a different format; 
support data are provided in the CCT-BT format (Section 2.1, 
Figure 2-2) and internal header records are provided in the 
CCT-AT and CCT-PT formats (Sections 2.2 and 2,3, Figure 2-12). 
Channels with defective data have been repaired in a manner that 
substitutes a minor image defect for a major one (Section 2.1.3). 
The cubic convolution resampling used to produce the CCT-PT 
images masks but does not eliminate the effects of duplicating 
lines. TM calibration on the intermediate CCT-AT product and the 
final CCT-PT product does not perfectly compensate for all 
problems, with the result that currently-available TM images show 
striping both by line and by scan (Section 2.2.5). Nominal 
geometric corrections which provide proper geodetic relationships 
between different parts of the image are available only for the 
final product image, the CCT-PT (Sections 2.3.2 and 2,3.3). 
Special user processing of the intermediate TM image formats can 
yield an image with approximate geometric alignment without 
resampling (Section 2.1.2); this special product, a CCT-B' ' or a 
CCT-A 1 ' image, is useful for radiometric studies, but does not 
otherwise take the place of the CCT-PT image. 


4.2 Research Areas 


During their investigation of TM image characteristics, the 
authors have noted many image or sensor characteristics that need 
to be determined and quantified. For this Early Results Sympo- 
sium, net all could be done; many interesting characteristics of 
TM images remain to be explored. 
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4.2.1 Overshoot & Undershoot 

The absence, or the presence and characterization, of 
undershoot and overshoot need to be determined for TM data. 
Kieffer et al. (1983) reported erratic overshoot and undershoot 
for a portion of the Washington D.C. scene (ID ■ 40109-15140). 
Since both cases represent the response lag of the sensors in 
crossing a sharp boundary in reflectivity, scenes featuring 
land-water or snow/ice-water boundaries need to be examined. 
Overshoot would exist if the scan, in passing from a uniform 
bright area to a uniform dark area, exhibits a gradational 
transition rather than a sharp transition. Undershoot would 
exist in the .case of a gradational transition from a uniform dark 
area to a uniform bright area. 

It is necessary to distinguish two cases. A gradational 
transition in the real world that is accurately recorded by TM 
sensors must be distinguished from a shafp transition in the real 
world that is inaccurately recorded by TM sensors as a grada- 
tional transition. For large areas, the boundary between bright 
and dark areas will be traversed by both forward and reverse TM 
scans. If the same gradational boundary is found in both scan 
directions, it should be a real-world event and not a sensor 
artifact, 

• Another consideration is the type of image product involved. 

Undershoot and overshoot on an image in B or A formats could be a 
characteristic of the sensor. The same effect on an image in P 
format could be a characteristic of the cubic convolution 
resampling procedure; this would be independent of forward and 
reverse scans, but should involve a narrow boundary zone. 


4.2.2 Bright-Spot Recovery 

Closely linked to the problem of overshoot and undershoot is 
the problem of bright-spot recovery. In the latter case, the 
item of concern is how quickly do TM detectors return to low 
reflectance conditions after sensing an individual pixel of high 
reflectivity . 

The evaluation done to date using specular reflectance 
pixels (Section 3.1) evaluated by printouts of DNs for small 
areas within an image indicates that the recovery from bright 
spots is sufficiently rapid. However, statistical studies using 
large numbers of specular reflectance pixels are needed to 
determine how fast the TM detectors recover. 


'4.2.3 Light Leaks 

The presence of internal reflections affecting the primary 
focal plane sensors was reported (Markham & Barker, 1984). 
Preliminary analysis indicates that the sensor receives light 
comming from two directions; from the front (towards the 
telescope) and from the back. Light that has passed through 
filters and sensors for band 1 has the greatest effect on other 
bands. In this case, there is a 15-pixel lag between sensing a 
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bright spot in band 1 and the effect of internally reflected 
light producing a 1.2% increase in DN value in band 4 (e.g,, a 
255 DN pixel in band 1 adds 3 DNs to what is sensed in band 4. 
The lag results from the different physical locations of the 
sensor array within the TM. There appear to be interactions 
among other bands, but at much lower percentages. 

The internal reflection will have to be completely 
characterized so that the effect can be eliminated before 
attempting to calculate the absolute radiometric data needed for 
a variety of research projects (Barker, Abrams, et al., 1984a). 


4.2.4 Out-of-Band Spectral Response 

Another possible phenomenon to be characterized is the • 
response of the sensors to light with a wavelength outside tne 
band width (Markham & Barker, 1983). In other words, the filters 
are not completely opaque outside of spectral window. Of concern 
is not only the light with wavelengths just outside the 50% 
shoulder that defines the spectral window for the filter, but 
also the light sensed that is at some wavelength distant from the 
filter window. How much light is sensed and at what wavelength 
location, both in absolute terms and relative to the filter 
window, are the questions for which answers need to be sought, 


4.2,5 Specifications 

How the TM system meets specifications (Anon., 1982b) will 
be discussed in a later document. This research is still 
ongoing. 


4.2.6 Odd-Even Difference 

Examination of image statistics has shown that there is an 
odd-even difference in detectors which is not removed by cali- 
bration. This results from a difference in physical location in 
the sensor array between even-numbered detectors and odd-numbered 
detectors (Anon., 1982b) and/or different electronics. This 
effect needs to be quantified and characterized using multiple 
images, so that it can be removed by improved radiometric cali- 
brations. 


4.2.7 Within-Scan Droop 

A possible problem is a change in the electronically-induced 
offset (Barker, Abrams, et al., 1984b; Metzler & Malila, 1983) 
during a scan. if the offset changes between the start and end 
of the individual scan, the usefulness of TM data for radiometric 
and digital-image classification-mapping purposes could be 
reduced. It should be possible to detect the effects of this 
condition, if it takes place, by the analysis of uniform bright 
and dark scenes. The effect, if present, would follow the 
direction of scan, and therefore could be separated from atmos- 
pheric and illumination effects. It should be possible to remove 
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the effects of this drift during radiometric calibration. 


* 


» 
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4.2.8 Spectral Calibration 

It is necessary to test the assumption of homogeneity of 
filters. The filters were cut out of a single piece of "uniform 
filter material (Markham & Barker, 1984), so that different 
filter chips came from different microlocations and small, 
detectable nonuniformities could be present. For characteriza- 
tion of the spectral calibration among data channels, it is 
necessary to look at the microscale heterogeneity of the filter 
material used. In addition, by comparing relative gains in each 
channel over several scenes, spectral differences between the 
channels can be characterized (Engel et al . , 1983; Markham & 
Barker, 1984). 


4.2.9 Autocorrelation 

The examination of the differences between nominal and IC 
calibration procedures has shown (Section 2.2.5) that changes in 
DNs between the two calibration procedures applied often occur 
in scan lines. Printouts of the difference images give ample 
visual evidence of autocorrelation (Figures 2-22 to 2-25). This 
needs to be tested statistically, and characterized so that auto 
correlation effects can be removed by those that consider it 
necessary to do so (Basu & Odell, 1974; Irons & Labovitz, 1982). 
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SECTION 5 — CONCLUSIONS 


5.1 Definitely True Conclusions 


Without any reservations or qualifications, we can make the 
following conclusions at this time about the Landsat-4 TM sensor 
system and its digital-data products: 

o Detector 3 of band 5 does not respond. In any image 
processing software system, this channel of data will 
need to be replaced {Section 2.1.3). 

o The "empty-bin" phenomenon in CCT-AT images results 
from integer truncations of mixed-mode arithmetic 
operations (Section 3.4). 


5.2 Probably True Conclusions 


With only minor reservations or qualifications, we can 
conclude : 

o TM final-product digital data (CCT-PT) is very useful 
data with high radiometric and geometric quality 
(Salomonson & Koffler, 1983; Bernstein et al . , 1983) 
that will open new avenues of Earth-related 'research 
(Williams, 1983 ) . 

o Interband registration, calculated from apparent 

spectral-response pixels for small, reflective objects 
on the ground, is within 0.5 pixel accuracy in both 
line and sample directions (Section 3.2), both within 
and between each focal plane. 

o Ancillary data files, residing on the same tapes as the 
image files, contain useful information in character 
format. This data can be printed in easy-to-read forms 
by utility "TAPEDUMP" programs that read and print 
character format. 
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5.3 Possibly True Conclusions 


With reservations and qualifications, we can make the 
following conclusions at this time: 

o The SCROUNGE-era fixup employed for the various 
defective data channels (Section 2.1.3), that of 
replacing the defective line with a duplicate of the 
previous good line, is acceptable. This conclusion 
should be tested by a study that compares accuracy 
against computer processing cost for these alternate 
methods: 1) complete line replacement using the 

previous line; 2) alternate pixel replacement from 
* previous and following lines; 3) replacement with 

the average of neighboring lines; and 4) calculation of 
missing scan lines from data in the same and 
neighboring lines in the different bands once band 
adjustments are made for differing spectral responses. 



I 


SECION 6 


RECOMMENDATIONS 


6.1 Landsat-4 Spacecraft Management 


As part of the continuing improvements and changes to the 
operational procedures for spacecraft control, maintenance, and 
sensor operation, we recommend that: 

o repeated lamp turn-ons be avoided when using 

manual-mode internal calibration in order to extend the 
life of the lamp (Barker, Abrams, et al . , 1984b). 


6.2 Landsat-4 SCROUNGE 


As part of the continuing improvements and changes to the 
SCROUNGE system, we recommend that: 

o Better documentation accompany the digital products. 

This includes not only documentation on paper but also 
internal documentation on magnetic tape files. More 
parameter values should be recorded in literal 
(character) fields so that a simple tape dump can 
provide the user with the data. 

o Better calibration be applied to remove channel and 
scan differences that result in striping (see also 
Bernstein & Lotspeich, 1983). This would involve 
treating forward and reverse scans separately rather 
than averaging pulse heights from forward and reverse 
scans (Barker, Abrams, et al., 1984b) and using 
drive-specific background average data rather than 
scene-wide averages to determine offset (Barker, Ball, 
et al. , 1984) . 

o Data from bad channels be handled better. For channel 
4 of band 2, the "bad" data should be used, not 
replaced. For channel 3. of band 5, a more appropriate 
substitution should be used (see also Bernstein & 
Lotspiech, 1983). 


6.3 Landsat-4 TIPS 


In the development of the TM Image Processing System (TIPS), 
we recommend that: 

o In addition to the SCROUNGE-era resampling by cubic 

convolution procedures to produce the CCT-PT image 
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(Section 2 , 3 . 2 ), the optional capability be developed 
to produce a CCT-PT product using the nearest neighbor 
resampling procedure. The latter procedure has the 
advantage of not changing the DNs of the pixels. Pixel 
values that are unchanged after tadiometric calibration 
are needed to allow radiometric characterization and 
evaluation of TM scenes, 

o In addition to the SCROUNGE-era geodetic correction to 
produce a CCT-PT product, TIPS apply corrections based 
upon known locations of ground-control points. 

Research into registration requirements for accuracy 
(Boyd et al., 1981) has shown that an even distribution 
of ground-control points is more important than the 
total number. 


6.4 Landsat-4 NOAA 


We recommend that: 

o Ground processing computer programs have the ability to 
go directly to CCT-PT geocorrected images. This would 
eliminate integer arithmetic problems with intermediate 
products, thus providing higher radiometric accuracy. 

o Images be made available to the user community both on 
a full scene basis and -on a partial scene basis. Such 
a partial scene basis could produce images that 
correspond to lx2-degree topographic quadrangles, or to 
15-minute quadrangles. Partial scene images could be 
made available on CCTs or on 8-inch floppy disks. 

o Archival, browse, and product images be recorded on 
.video laser disk media. Low cost, random access, 
•multiple scene storage, and ease of distribution are 
the obvious advantages (Goldstein, 1982; Mobery & 
Laefsky, 1982; Moussy, 1983). Remote sensing has been 
described as a "write-once , read often" process. 


6.4 Landsat-D' 


In preparing the TM sensor system before the launch of 
Landsat-D', we recommend that: 

o The dynamic range be expanded. Alternately, high-gain 
/ normal-gain options should be added to the 
analog-digital converter to permit a normal scene to 
more nearly fill the available DN range of 0-255. TM 
bands 2 and 3 do not reach saturation levels, even when 
recording light reflected from snow. 
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o 


An opaque light-trap be placed behind the sensors to 
eliminate internally reflected light which is affecting 
the TM sensors of Landsat-4. 


o An opaque light-mask be placed around filter mounts on 
the prime focal plane to eliminate light passing 
between the detector assemblies (Barker, Abrams, et 
al., 1984b) . 
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Figure 2 4. An example of a Band 5 image in CCT— AT format. The "before" image demon 
strates the existing scan shift between forward and reverse scans (See Figure 2—5). 
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Figure 2-7. An example of data for TM band 5, showing the defective channel 3. The scan 
organization is the same shown in Figure 2-3. 


r 


ORIGINAL 



i 

I 


O' 

c 

O I 
— t r\ 
3 

<0 
V ) v 
H a> u 
a E 3 

I "H (J\ 

f-* 4J ■-* 

U U. 

u 

•v .C 

"O W *J 

0) *H 

.x u 5 

O 3 

ooo) 

— < U-l U 

13 n 
C T 3 0. 
3 O E 
4J O 
c 3 u 

a 

* Q • 

'V U O 


O' 


*T3 0) 

3 


C U 

E 


3 <y 

• 


n 5 

4J 


s u 

3 


E- 0 

a 


10 

u 


u c 

0 


0 0 4-1 


4-1 (0 

a 


C O 

o -c 

o 

• 

■H 4J 

a 


4-> 

3 E 

4J 


N 0 

0) 


•H u 

u 


C 4~i 

3 


3 

TD 


O' to 

0 


Ij ~l 

u 


0 0) 

a 


X 



T3 

0 


C Oil 


3 



> 

O' 


3 O 

c 


4-> Cid 

■H 


fTJ >-l 

10 


t3 

to 


E 

o 


U-l o 

u 


0 -N 

o 


H 

u 


0) 

a 


3 

0) 


a 3 

u 


E u 

a 


3 

X 4-1 

73 


0) 3 

Q 


J 3 

a 


C 4J 
< 

< 


o 

o 


4-1 



• o 

4J 


a 2 

l 

>1 

• 

,~N| 

£3 

■D 

• 


*— 1 

0) 4-J 

O 

1 

Ij 3 

C 

JN 

3 E 

r~4 


O' U 

H 

X) 

•H 0 


c 

Cu 4-1 

3 

3 



















LANDSAT-4 TM - DIGITAL TAPE CONTENTS - CCT-AT 

(4X BLOCKING) 

SoROUNGE-ERA 11982-1983) 


oricin;.;. r> ;; 

™ POOR Q’J ALllY 


V) 

3 

>1 

X) 

T3 

<D 

W 

3 





• a 

0 u 

■ — < r} 

1 

rs nj 
■u 
O <T3 

3 

O' U 

o 

X *4-t 




LANDSAT-4 TM CCT-AT IMAGE FILE CONTENTS 
SCROUNGE ERA (1982-1983) 


ORIGINAL PAt»S 10 
OF POOR QUALITY 


v 


fi 



LANDSAT-4 TM - CCT-AT - INITIAL PORTION OF DATA 

RAATRSX (UNBLOCKED) 

MEMPHIS TN (ID = 40037-16033) 22 AUG 1982 


ORIGINAL PAQii B 

OF POO'? miAUTY 


U3awnN 13NNVHS 








ct 

rt 

* 

« 

PI 


o 

O' 

a 

fv 

CT 

CT 

* 

« 

PI 

mm 







mm 

mm 

mm 


*“» 













up 

a 

pi 

PI 

PI 

O 

N 

a 

fv 

o 

« 

N 

a 

CT 

O 

CT 

« 

CT 

sO 

o 

pi 


(M 

o 

« 

CT 

CT 

w 

W 

n 

CT 

OS 

© 

a 

a 

© 

0*4 

© 

O 

O 

© 

o 

o 


* 

IP 

Cl 

CJ 

PI 

N 

© 

N 

CK 

© 

© 


CO 

o 

o 

CP 

« 

,o 

O' 

(It 

Vl 


N 

* 

CT 

CT 

CT 

CT 

so 

w 

CT 

OS 

C0 

o 

CO 

© 

© 

o 

o 

co 

CO 

o 

© 


o 

* 

d 

CJ 

PI 

UP 

© 

pi 

CO 

cs 

(S 

to 

© 


o 

o 

© 

© 

© 

© 

o 


PI 

co 

n 

CT 

CT 

w 

so 

SO 

SO 

© 

o> 

CO 

© 

CO 

© 

a 

to 

CO 

CO 

CO 

© 


Cl 


Cl 

PI 

PI 

IP 

© 

N 

N 

•* 

N 

O' 

pi 


© 


o 

© 

© 

Q 

© 


N 

£ 

n 

CT 

« 

CT 

SO 

SO 

CT 

CO 

CO 

CO 

O' 

o 

© 

© 

© 

© 

a* 

© 

© 



N 

Cl 

CJ 

Cl 


SO 


N 

CT 

CO 

co 

© 

•M 

N 

a 


CO 

CO 

o 

PI 


Cl 

so 

cp 

CT 

CT 

CT 

CT 

N 

*o 

N 

CO 

5 

© 

© 

CO 

CO 

co 

CO 

a 

© 

© 


© 

UP 

Cl 

PI 

PI 

Cl 

SO 

N 

UP 

O 

O' 

CO 

© 

SO 

S' 


PI 

so 

CT 

o 

fv 


PI 

* 

w 

CT 

« 

so 

SO 

CT 

CT 

N 

CO 

a 

© 

CO 

o 

o 

o 

© 

© 

a 

co 


© 

o 

PI 

PI 

PI 



W 

CT 

CO 

PI 

2 

© 

vO 

SO 

■co 

O' 

CT 

VO 


© 



NO 

CP 

O 

CT 

CT 

SO 

CT 

SO 

SO 

O' 

a 

a 

O 

a 

a 

© 

o 

© 

CO 

CS 


c 

N 

Cl 

PI 

Pi 

W 

»0 

SO 

N 

o 

UP 

© 

sO 

© 

Cl 

pi 

CP 

« 

CP 


© 



O 

cp 

CT 

« 

CT 

SO 

sO 

SO 

tv 

CO 

CO 

CO 

fv 

a 

I? 

CO 

CO 

CO 

a 

CO 


N 

N 

pi 

CJ 

PI 


?i 

N 


* 

* 

N 

eo 

PI 

© 

© 

2 

n 

CP 

a 

o 


mm 

VO 

cp 

CT 

CT 

t» 

tv 

CT 

tv 

N 

tv 

N 

tv 

a 

a 

CO 

G3 

a 

G3 

N 

N 























00 

CO 

SO 

Cl 

<N 

PI 

PI 

*■* 

a 

PI 

N 

N 

w 

W 

sO 

© 

© 


Pi 

o 

© 

© 

© 

Cl 


cp 

CT 

CT 

CT 

N 

SO 

P- 

tv 

N 

N 

N 

N 

CO 

CO 

CO 

CO 

a 

CO 

t- 

cs 

, T-* 






















cli 

CO 

up 

M 

CJ 

PI 

PI 

o 

N 

N 

CT 

M 

PI 

CT 

© 

© 


o- 

© 

pi 

PI 

CO 

© 


n 

CT 

CT 

CT 

N 

sO 

tv 

fv 

N 

tv 

N 

N 

N 

a 

a 

a 

£3 

ca 

N 

a 

07 


C4 

Cl 

PI 

CJ 

N 

o> 

PI 

* 

SO 

PI 

W 

© 

** 

o 

*~4 

© 


CTPf 

CO 

o 



CP 

CT 

CT 

CT 

O 

SO 

N 

N 

N 

tv 

fv 

fv 

CP 

s 

a 

CO 

o 

CO 

N 

cs 

< 

CP 

IP 

CJ 

CJ 

PI 

© 

CO 

SO 

O 

"C 

IS 

N 

© 

CT 

CP 

a 

o 

sO 

SO 

N 

CP 

K 


o 

CT 

ra 

CT 

so 

SO 

fv 

N 

N 

N 

t» 

CO 

CO 

CO 

s 

s 

N 

N 

tv 

cs 

LU 

04 

© 

CJ 

CJ 

CJ 

M 


sO 

Os 

W 

'P 

© 

f 

UP 

CP 

a 

CO 

O 

e 

© 


Ui 



CT 

CP 

CT 

tv 

tv 

N 

N 

N 

t- 

CO 

a 

CO 

CO 

N 

s 

CO 

CO 

o 

cs 

iD 



CJ 

Cl 

CJ 

bp 


up 

OS 

* 


<#■ 

n 

-c 

s* 

O' 

a 

a 

CO 

CT 

* 

2! 

M 

6 

CP 

CP 

CT 

N 

N 

N 

t- 

tv 

CO 

a 

a 

a 

CO 

a 

s 

t* 

N 

CO 

cs 

, D 


•"* 












• 









o 

o 

CJ 

Cl 

CJ 

* 

CJ 

CT 

© 

© 

SO 

■Sf* 

* 

SO 

w 

© 

PI 

© 

o 


IP 

O 



CP 

CT 

CP 

N 

N 

N 

C3 

CP 

CO 

CO 

a 

a 

CO 

CO 

CO 

a 

CO 

CO 

cs 

cs 






















o 

© 

© 

Cl 

c» 

CJ 

rC 

Q 

CP 

CS 

Pi 


w 

UP 


UP 

UP 

CO 

CP 

CP 

fv 

t- 

w 



CP 

• 19 

CT 

N 

N 

tv 

fv 

a 

CO 

C3 

a 

CO 

a 

a 

s 

Q 

a 

cs 

cs 

C5 

C3 

CJ 

CJ 

Cl 

UP 

O 

© 

O' 

© 

* 


* 

UP 

N 

UP 

CO 

CO 

CO 


Pi 




CT 

CT 

CT 

t- 

t- 

vO 

N 


a 

Q 

a 

CO 

a 

CO 

s 

N 

N 

cs 

CS 



o 

CJ 

Cl 

CJ 

UP 



CP 

© 

© 

UP 

O' 

CT 

N 


CT 

CO 

CO 

CT 

O' 

' 



CP 

« 

CT 

N 

N 

t- 

a 

© 

mm 

© 

CO 

a 

CO 

a 

a 

N 

N 

t* 

N 

cs 


vO 

Cl 

CJ 

CJ 

CJ 

N 

CJ 

* 

CP 

N 



SO 

© 

SO 

SO 

SO 

VO 

sO 

t- 

PJ 



© 

CT 

CT 

CT 

t* 

N 

N 

CO 

© 

© 

P5 

a 

a 

a 

CO 

CO 

c z 

CO 

N 

cs 


Up 

c? 

Cl 

CJ 

CJ 

N 


Os 

Os 

© 

w 

UP 

pi 

w 

w 

IP 


fmO 


PI 

* 



so 

CP 

CP 

CT 

N 

t- 

N 

CO 

© 

© 

a 

© 

© 

© 

© 

o 

o 

Q 

o 

© 

















•■n 



** 

** 





CJ 

Cl 

CJ 

S3 

N 

a 

a 

© 

N 

a 

P3 

M 

CO 

N 

Sf 

o 

CP 

© 

CO 




CT 

CT 

CT 

CO 

a 

a 

a 

© 

*■* 

CO 

CO 

© 

o 

z 

PI 

CT* 

c 

O 

© 

CS 


CP 

© 

CJ 

CJ 

CJ 

SO 

SC 

pi 

sO 

N 

w 

CP 

O 

CT 

© 

UP 

SO 

© 

o 

CO 

IP 




CT 

w 

CT 

a 

a 

a 

CO 

CO 

a 

© 

o 

© 

© 

© 

CO 

C3 

a 

N 

N 


CJ 

o 

Cl 

CJ 

CJ 

r* 

UP 

Os 

up 

•o 

tv 

tv 

o 

CP 

t- 

VO 

S 

N 

N 

mm 

© 




CT 

n 

CT 

a 

CO 

a 

© 

mm 

CO 

CO 

O' 

a 

CO 

s 

tv 

S 

N 

N 

N 

N 



© 

CJ 

Cl 

CJ 

SO 


» 

© 

© 

CT 

© 

sO 

sO 

CP 

PI 


PI 

pi 

© 

_ 




CT 

CT 

CP 

0> 

O' 

a 

© 

N 

tv 

tv 

fv 

t» 

N 

S 

s 

N 

N 


h* 


w 

as 

- 

CM 

CT 


w 

CT 

N 

CO 

© 

© 

I* 

PI 

CT 

f 

CT 

so 

h* 

CO 

•*« 

© 

20 


Li 






















w 

j 

U3QV3H 






Viva 3DV1AII 








1 




1 


















w 

W 0) 

ra p 

H 

** (0 • 


Figure 2-12. An example of data and ogranization for TM images in CCT-AT format 
unblocked when written to a copy tape or a disk file. The three lines of "32" v 
represent blanks padding the internal header record before the data is unblocked 
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figure 2 13. An example of TM data in CCT— AT format, unblocked. Note the changes in DNs 
between forward and reverse scans. Compare with Figures 2-14 to 2-16 to see the same 
pattern in other multispectral bands for this same image segment. 
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Figure 2—14. An example of TM data for spectral band 2 in CCT— AT format, unblocked. Note 
the duplication of DNs for channel 5 to replace defective data in channel 4. Note also 
the difference in DNs between forward and reverse scans resulting from the scan offset. 
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Figure 2-15. An example of TM data for spectral band 5 in CCT-AT format, unblocked. Note 
the duplication of DNs for channel 4 to replace defective data in channel 3. Note also 
the difference in DNs between forward and reverse scans resulting from the scan offset. 
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Figure 2-18. An example of TM data for spectral band 2 in CCT-AT format, unblocked, 
internal calibration. Note the duplication of D*Ns for channel 5 to replace defective data 
in channel 4. Note also the difference in DNs between forward and reverse scans resulting 
from the scan offset. Compare with Figures 2-14 and 2-22. 
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Figure 2-20. An example of TM data for spectral band 5 in CCT-AT format, unblocked, 
internal calibration. Note the duplication of. DNs for channel 4 to replace defective data 
in channel 3. Note also the difference in DNs between forward and reverse scans resulting 
from the scan offset. Compare with Figures 2-15 and 2-24. 
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Figure 2 22. An example of the difference image between TM data processed using two 
different calibration options. Compare with Figures 2-14 and 2-18. 
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Figure 2-23. An example of the difference image between TM data processed using two 
different calibration options. Compare with Figures 2-13 and 2-23. 
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Figure 2-24. An example of the difference image between Trt data processed using two 
different calibration options. Compare with Figures 2-15 and 2-24. 
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Figure 2-25. An example ot the difference image between TM data processed using two 
different calibration options. 


ORIGINAL ! 

OF POOR Qi 


• • 


BAND AVERAGE COMPARISON OF DIFFERENT 
CALIBRATION METHODS 

LANDSAT-4 TM - CCT-AT - NOMINAL MINUS 1C HISTOGRAMS 
MEMPHIS TN (ID - 40037-160331 22 AUG 1982 
START LINE (SL) - 3001, START SAMPLE (SS) - 2001 NL - 1024 

NS - 1024 

SCROUNGE ERA (1982-1983) 


LEAST 

VALUE 

FRACTION 

-3 

.000001 

-2 

.000564 

-1 

.331877 

0 

.666902 

1 

.000657 

-3 

.000001 

-2 

.000025 

-1 

.140252 

0 

.845494 

1 

.014228 

-2 

.000001 

-1 

.159622 

0 

.834955 

1 

.005421 

2 

. 000**7 

0 

.083167 

1 

.864659 

2 

.052174 

1 

.037438 

0 

.912251 

1 

.050310 

0 

1.000000 

-2 

.000899 

-1 

.231982 

0 

.702545 

1 

.063312 

2 

.001261 

3 

.000000 

4 

.000000 

S 

.000000 

6 

.000001 


HISTOGRAM 




Figure 2-26. Histograms showing the frequency of differences 
between images produced by two different TM calibration options. 
Compare with Figures 2-22 to 2-25. 
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Figure 2-27. Histograms for each channel showing the different frequency of differences 
between band-7 images produced by two different calibration options. Compare with Figure 
2-2b . 
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Figure 3-4. An example of a histogram for one data channel showing empty bins at DN 
values of 18, 51, 84, & 117. 
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example of empty bin locations and periodicity for data from detector 


